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ABSTRACT 


Methods  are  well  known  for  determining  testing  times  to  minimize  the 
mean  cost  of  testing  plus  mean  cost  of  an  undetected  failure  (linear  in  the 
mean  time  between  failure  and  detection),  when  testing  does  not  degrade  a 
good  system.  Here,  we  introduce  a model  in  which  the  ith  test  may  either 
cause  a failure,  with  probability  j3 , or  increase  the  remaining  failure  rate  to 
*1  >^i_1  without  changing  the  form  of  the  conditional  lifetime  distribution. 
Algorithms  are  given  for  finding  the  best  testing  times  in  cases  of  uniform 
and  exponential  failure  time  distributions.  Optimization  over  a single 
cycle  is  considered  first,  and  then  the  case  with  component  renewals  is 
solved  using  the  mean  loss  per  unit  time  criterion. 
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1.  Introduction 

Barlow,  Hunter  and  Proschan  ^ posed  and  solved  an  interesting 

class  of  inspection  problems  related  to  reliability  and  maintainability  of 
systems.  In  their  models,  a system  operates  for  random  time  until  it  fails, 
but  the  failure  can  be  detected  only  after  a test  costing  Cj  has  been  per- 
formed. For  example,  judgment  of  whether  or  not  a radar  system’3  de- 
tection and  false-alarm  probabilities  are  acceptable  might  require  connec- 
tion of  special  test  equipment  and  injection  of  special  testing  signals. 

One  approach  to  this  problem  seeks  optimal  testing  times  to  minimize 
a loss  function  of  the  form 

L = E jcj  N + c2  dj  (1.1) 

in  which  N is  the  number  of  tests  until  the  first  one  after  the  failure,  and  d 
is  the  time  between  the  failure  time  and  the  subsequent  detection  time  t^  . 
Reference  3 gives  algorithms  for  finding  the  optimal  inspection  times  t^  for 
lifetime  distributions  which  are  uniform,  exponential  or  Polya  of  order  2, 
when  using  the  single-cycle  (no  renewal)  loss  function  in  (1.  1). 

A second  point  of  view  postulates  repair  or  replacement  after  each 
failure  detection,  with  associated  average  cost  and  average  repair  time  re- 
quirements. In  this  case  the  performance  is  judged  using  the  mean  cost  per 

(14) 

unit  time  or  equivalently,  the  ratio  of  mean-cost-pe r- renewal  to  mean- 
time-between- renewals.  Reference  3 shows  how  techniques  used  for  the 
single-cycle  problem  can  be  modified  to  find  the  best  testing  times  for  the 
case  with  renewals. 

Here  we  seek  similar  solutions  but  include  the  possibility  that  the  me- 
chanical and  electrical  stresses  of  the  test  might  destory  or  damage  the 
system.  Mathematically,  we  allow  the  ith  test  to  destroy  the  system  with 


-3- 


probability  0 or,  with  probability  (1-0),  to  increase  the  failure  rate  (re- 
ciprocal of  mean  time-to-failure ) to  X^  > X^  1 changing  the  form 

of  the  conditional  lifetime  distribution.  These  generalizations  are  approached 
in  a direct  manner  by  using  dynamic  programming  to  set  up  recursive  opti- 
mization equations.  The  main  contribution  here  is  the  presentation  of  con- 
vergent algorithms  for  solving  the  optimization  equations. 

Section  2 examines  exponential  conditional  lifetime  distributions,  and 
Section  3 gives  similar  results  for  uniform  conditional  distributions.  Sec- 
tion 4 discusses  extensions  to  more  general  distributions,  etc. 


2,  Exponential  Conditional  Lifetime  Distributions 

We  consider  a single  component  system  with  two  conditions  good  or 
failed.  The  system's  status  can  be  determined  only  by  a test  costing  c^ 
units  a failed  system  introduces  a penalty  of  c ^ per  unit  time  during  the  in- 
terval between  the  failure  time  t and  the  subsequent  inspection  time  t^. 

The  inspection  times  (to  be  optimized)  are  denoted  by  t^,  t^,..,,.  with  in- 
terchecking times  defined  as  6 = (t  , , - t ). 

fC  /Ct  * K. 

The  exponential  conditional  lifetime  assumption  enters  through  the 
family  of  conditional  failure  time  densities 

. -X  (t-t  ) 

(i)  f (t]t  <.t)  = X e K K ; t >t 
K K K 


mi  <t<tJt+I>. 


X e 

K 


-X  (t-t  ) 

K ' K' 


(2.1) 


. - X 6 

1-e  k K 


While  any  sequence  of  test-dependent  failure  rate  parameters  might  be  used, 
the  algorithms  to  be  used  below  require  that  X — « oo  as  the  number  k of 
possible  tests  increases  without  bound. 
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2.  1 Single-Cycle  Exponential  Case 


We  consider  first  the  case  in  which  a test  never  causes  immediate 
failure  ((3=0)  and  in  which  performance  is  judged  over  a single  failure  cycle. 
When  including  test  induced  performance  degradation  it  is  reasonable  to 
add  another  cost  function  reward  term  (-c  per  unit  time  of  good  operation. 
Otherwise  an  optimal  policy  might  call  for  several  quick  tests  which  in- 
crease the  probability  of  short  lifetime  but  also  decrease  the  mean  time  be- 
tween failure  and  detection. 

With  N as  the  random  variable  denoting  the  index  of  the  test  which 
actually  detects  a failure,  the  loss  to  be  minimized  by  the  choice  of  testing 
times  can  be  expressed  as 

L0(trt2...  ) = e|Ci  N + c2(tN"  t)  - c3(t_t0)  I t>t:a]  (2<2) 


Since  the  testing  times  after  t can  influence  only  future  contribution  to  the 
mean  cost,  it  is  useful  to  consider  the  optimal  mean  future  loss  L,°  at  any 
testing  time  t^.  This  can  be  expressed  recursively  as 


L ° = min  (1-e  K 
* t 

K+l 


\ 6.  A+i 


K)  f [C1  + c2  (tK+l  ■ fc)  ■ c3(t"tK)]  £ K 


-\  6 
K K 


[C1  + c2  ( Sc+  1 " fc)  + LK  + 1. 


(2.  3) 


in  terms  of  the  two  cases  t <t<t  or  t > t . , . 

K K+l  K+l 


The  expression  for  in  (2.  3)  simiplifies  to 


L°  = min  Cj  - (c2+c3)Ak+  c2  + e 
0 

K 


(2.  4) 
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A recursive  equation  for  the  optimal  interchecking  time  6°  follows  from 
simple  differentiation. 


,0  — 1 g /C  T O . I 

5/c  " \ ^ <+  i + 1 + C3/c2 


(2.  5) 


Substitution  of  (2.  5)  into  (2.4)  yeilds  the  companion  relation 


L°K  = C1  - c3/Xk  + c2 


(2.6) 


At  first  glance  (2.  5)  may  not  appear  to  be  useful  necessary  condition, 

since  it  is  defined  in  terms  of  the  next  optimal  loss,  which  is  as  yet  unknown. 

Futhermore,  the  structure  of  the  problem  does  not  permit  definition  of  a 

"last"  step  from  which  we  can  work  backwards.  However,  the  assumed 

growth  of  X^  does  permit  indirect  evaluation  of  the  "last  step"  optimal  cost 

as  c j,  as  follows . As  k and  X ^ grow  very  large,  P[  t<t  +A  [t>t  ] — • 1 

for  any  A>  0.  Thus,  in  that  limit,  t°  . . —»t  and  L°  c . 

More  formally, 

Lemma  1.  lim  L°  = c . 

K — oo  k 1 

Proof;  The  future  mean  loss  for  any  (not  necessarily  optimal)  sequence  of 
testing  times  may  be  written  similarly  to  (2.4)  as 

-X  6 

Ln(5fc'  *•*  } = cl"  (c2+c3,/X/c  + c2  6/c+e  K 'C(LK+1+(c2+C3)/XK)  (2‘7) 


The  limiting  expression 


K he  “ c , + e2  5, 


is  clearly  minimized  by  5^  = 0,  with  the  corresponding  limiting  loss  stated 
in  the  lemma. 
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The  foregoing  boundary  condition  on  L°  suggests  the  following  algorithm 
for  computation  of  the  optimal  testing  intervals  in  relation  to  a given  failure 

rate  sequence  X q 

1.  Choose  an  initial  maximum  number  of  tests,  Mj  . 

2.  Assume  L°(Mj)  = c j and  use  (2.  5)  and  (2.  6)  to  compute 


jlMj), 


LM1-1(M1)* 


-2^1^’ ' * ’ ' 6 0^M1^’ 


L°(Mi) 


3.  Repeat  for  Mg  > M^,  etc. 

4.  Stop  when  6°  (M<+1)  and  6°  (MJ  are  sufficiently  close  in  respective 
values  for  1 < i <_I,  and  a suitable  I.  (I  is  specified  with  respect  to  a 

sufficiently  high  value  for  P[  t < tj.  ] ) 

With  regard  to  sufficiency  of  the  zero-derivative  condition  underlying 

(2.  5),  it  is  easy  to  show  that  6^  > 0 implies 


(2.8) 


and  that  (2.  9)  insures  that  the  second  derivative  of  the  bracketed  term  in 


(2.  4)  is  positive. 


Examples 

1.  The  first  example  has  c^  = 1,  Cg  = 20,  c3  = 20,  \q  = 2 and 
X,c  = X0/(0.9)’< 

i. e each  test  reduces  the  mean  residual  lifetime  by  a factor  of  0.  9.  Table  1 
shows  the  approximate  values  for  6°  and  L°  which  result  from  - 21. 
Similar  calculations  with  Mg  = 31,  M3  = 41  etc  produced  the  same  10  digits 
for  the  first  seven  6.  and  L.  values.  This  degree  of  convergence  will  be 
adequate  if  the  resulting  P [ t < t?]  * 1 - exp  [ - g 6.  X . ] = 0.  98  + is  ac- 
ceptable. 
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Table  1 also  shows  how  the  tests  reduce  the  mean  lifetime  from  0.  5 
(in  the  absence  of  testing)  to  0.  438.  The  EJT)  shown  there  is  defined  as 
the  mean  lifetime  given  exactly  k tests,  using  the  6^  values  found  above. 

This  assumes,  for  the  purpose  of  approximate  caluclations,  that  a failure  after 
time  tK  goes  undetected.  The  convergent  E^  sequence  can  be  computed  re- 
cursively from  the  relation  ^ 


E|c,t)  = Ejc_i(T)-  ( f 

X-  1 


K-  1 

, - 2 \ . 5. 

T-  )e  0 11 

K 


Similar  calculations  with  other  parameter  values  reveal  interesting 
trends.  Interchecking  times  are  smaller  if  \ is  increased,  if  c„  is  in- 

w u 

creased,  or  if  c3  is  decreased.  If  c3  is  large  compared  to  c and  c then 
there  is  a premium  on  long  life  and  L°  is  a monotonically  increasing  func- 
tion of  initial  failure  rate.  This  is  demonstrated  in  Figure  1 where  it  should 
be  noted  that  , e.g.  , the  loss  after  tj,  when  having  started  with  XQ  at  tQ,  is 
equivalent  to  the  loss  after  tQ  when  starting  with  (XQ/p).  (L^Xg)  = LQ(\0/p)) 
However,  Figure  2 shows  that  if  c 3 is  relatively  smaller,  then  increasing 
the  initial  failure  rate  may  increase  or  decrease  the  mean  loss. 

Table  2 shows  results  analogous  to  those  in  Table  1 except  that  here 
the  failure  rate  grows  linearly  instead  of  geometrically: 

= X0  ( 1+K  > 

— Optimum  Inspection  with  Renewals  - Exponential  Cas 
The  foregoing  model  will  now  be  generalized  to  allow  for  repairs  or 
renewals  which  each  have  a mean  cost  s and  a mean  completion  time  of  r 
time  units.  Inspections  times  between  renewals  must  again  be  selected, 
but  performance  is  now  measured  in  terms  of  the  ave 


rage  cost  per  unit 


time.  The  cycle,  starting  at  time  tg  and  ending  at  t^j  , will  accrue  a cost 

i 1 

C.  = N.  c.  + d.  c„  + s (2.  9) 

l l 1 l 2 

and  will  have  a duration 


T . = (h.  - t ) + r 
l N.  0- 

l i 


(2.  10) 


Using  these  terms  the  average  cost  is  defined  as 


k 

S C. 

l 

R(6)  = Lim  - * 


(2.  11) 


The  argument  6 is  shown  to  emphasize  the  dependence  of  this  average  cost 

on  an  inspections  policy  6=  j 6^,  6 ^ ...  V . 

(4) 

It  is  well  known  that  R(6)  can  be  expressed  in  terms  of  the  mean 
cost  per  cycle  C and  the  mean  cycle  duration  t,  as 


R( 6)  = C / t (2.12) 

Both  the  numerator  and  demonenator  of  (2.  12)  are  affected  by  the  inspection 
policy  6.  However,  it  is  possible  to  find  the  6 which  minimizes  R(6)  by  con- 
sidering an  auxiliary  optimization  problem  having  the  cost  function 

£ (p,  6)  - c - u T (2.  13) 


These  two  problems  are  related  as  follows. 

(3)  * 

Theorem  1:  If  there  exists  a |j  = u for  which 


min 

6 

imizes  R(5). 


# A 

JC  (u  ,6)  = 


£ (u\6(  n ) ) 


0,  then  the  schedule  6 ( y ) also  min- 


-9- 


The  existance  of  such  a desired  p*  is  quite  evident  for  the  present 
problem.  Equation  (2.  13)  can  be  rewritten  as 


L ( y,  5)  = c N + c^  d+s-^[  T + r] 


(2.  H) 


a little  thought  shows  that  when  p _>  c7  , increases  in  6.  will  decrease  £ 
since  t > d.  Thus,  for  p>  c^  the  optimal  policy  is  to  make  no  inspections 
or  repairs,  with  a corresponding  optimal  £ < 0.  On  the  other  hand,  when 

P = 0,  £ >0  for  any  policy including  the  optimal  one.  We  conclude  that 

£ (p,  6, (p))  changes  from  a positive  value  at  p = 0 to  a negative  one  for  p = c7, 
and  this  change  must  be  a continuous  one  passing  through  zero  at  least 
once  in  view  of  the  smooth  nature  of  all  densities  and  cost  terms  entering 
into  the  loss  functions. 

An  algorithm  which  uses  Theorem  1 to  optimize  inspection  time  in  the 
presence  of  renewals  consists  of  two  parts.  An  algorithm  like  the  one  in 
Section  2.  1 is  needed  to  minimize  (2.  14)  for  each  choice  of  p until  £ (4,  6( p) ) 

= 0.  The  first  part  parallels  the  approach  in  the  single-cycle  case  discussed 
earlier  since  £ defined  in  (2.  14)  is  a single  cycle  loss  with  p similar  to  c^. 

In  particular  (2.14)  leads  to  the  following  recursive  expression  for  the 
minimal  mean  future  loss  after  test  jc,  ^ 


L°  = min  jcj  + (c,  - p ) 6^  - c2/\k  + ( s - u r) 

I 


6 

K.  K.  T O 


Differentiation  yields  the  recursive  optimization  equations. 


-10- 


(i) 


c ,+X 

2 K 


'k+1 


\(s  - 4 r)  ) / (c2  - p) 


(2.  16) 


(ii)  = ci  “ dAK  + ( s " H r)  + (c2  - p)  6° 

It  can  be  shown  ^ that  the  p < c2  condition  derived  above  insures  that  the 
zero- derivative  condition  corresponds  to  a minimum  and  that  all  6°  > 0. 

Finally,  Lemma  1 must  be  modified  for  the  present  loss  function  to 

Lemma  2:  ^lirn.  £ ° = c ^ - p r + s (2.  17) 

Once  the  p of  Theorem  1 has  been  found,  it  follows  immediately  that 


min  R ( 6 ) = 
6 


* 

d 


(2.  18) 


Sensitivity 

The  preceding  paragraphs  describe  how  to  get  the  best  inspection 
schedule  6 and  the  corresponding  minimal  loss  per  unit  time.  It  is  also  of 
interest  to  see  how  sensitive  this  solution  procedure  is  to  the  precision  of 
knowledge  of  model  parameters,  e.g.  Xq,  and  how  sensitive  it  is  to  the 
accuracy  of  the  iterative  calculations . It  is  straight  forward  to  use  the  ex- 
ponential conditional  densities  to  get  interative  expressions  for  N,  d and  t 
which  can  be  terminated  when  the  additional  terms  contribute  very  little. 

(See  for  example,  the  expression  in  the  example  of  Section  2.  1)^ 
Examples  Tables  3,4  and  5 show  results  for  numerical  examples  based  on 
the  models  in  this  section.  Table  3 shows  the  small  sensitivity  of  R(6)  to  use 

i * 

of  non-optimal  testing  schedules  that  result  when  p d • The  parameters 
for  Tables  4 and  5 differ  only  in  that  the  latter  has  non- zero  values  for  r 
and  s,  the  renewal  costs.  Both  of  those  tables  show  minimum  R values  for 
severed  X q as  well  as  the  performance  degradation  when  X ^ = 5 but,  when  the 
optimum  6-schedule  for  a different  value  is  used. 
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3.  Uniform  Conditional  Distributions 

In  this  section  we  consider  a failure  model  in  which  a single-component 
system  begins  operation  at  tg  with  a uniform  failure  time  distribution  between 
tQ  and  (tg  + T0).  An  inspection  test  at  time  tj  = (tQ  + 6j)  of  a "good"  system 
leaves  a remaining  failure  distribution  uniform  between  t^  and  t^  + where 
Tj  = p (Tq  - 6j),  G<  p < 1.  That  is,  after  each  test  the  conditional  distri- 
bution is  still  uniform,  but  with  its  mean  lifetime  reduced  by  a factor  p,  as 
compared  to  the  case  of  no  test. 

The  analysis  of  this  uniform  conditional  distribution  model  is  not  so 
neat  as  that  of  exponential  models,  but  convergent  algorithms  have  been  found 
for  determining  the  optimum  number  and  locations  of  inspection  tests  for 
given  to  and  p when  the  single-cycle  loss  function  of  (Z.2)  is  to  be  minimized. 

Dynamic  programming  is  used  to  optimize  the  location  of  the  first  of  n 
tests  (tj  (n)  = tg  + 6j(n))  assuming  that  all  future  tests  are  optimally  located. 
Difficulties  arise  when  a zero-derivative  condition  yields  an  inadmissible 
5j(n)>  Tq  or  when  the  remaining  lifetime  T^  is  so  short  that  the  analytical 
minimal  future  loss  expression  L°n  j (T  ^ ) u^edfor  computing  6j(n)  is  no  longer 
valid.  These  pathologies  and  techniques  for  circumventing  them  are  best  des- 
cribed by  the  following  small  example. 

The  best  location  for  a simgle  test  is  found  by  minimizing  the  following 
mean  loss  expression  with  respect  to  x^  ( x°  = 6j).  (We  introduce  the  extra 
symbol  x to  help  distinguish  between  zero-derivative  intervals  and  optimal  in- 
tervals ). 

*1 

L1  (T0’  xl)  = C1  + T^  fQ  xj  l cl(xl  ‘ ‘ c3t  J dt 

x 1 <T0  - x . ) 

+ ( 1 - f~  ) [ - c3  Xl  + p ? (c,  - c3  ) ] (3.  1) 
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This  reduces  to 


L1(T0’  xl)  = X1  (c2  + c3  “ a0}/2  T0  " Xl(c3  " a0)  + C1  * a0  To/2  (3,2) 


in  which  we  use 


a0=  P (c3-  c2) 


(3.  3) 


The  minimum  of  the  parabolic  function  (3.2)  is  clearly  achieved  by  x^(l) 


c , - a 

x.  (1)  = 2_  T _ § (i) 

1 C2+C3"a0  0 1 


(3.4) 


which  lies  between  0 and  TQ  since  c3  - aQ  > 0.  The  corresponding 
mum  loss  is 


mim- 


T o /n,  , _ a0  T0  T0(c3-a0) 

L1  (T0>  “ c l 2 " 


2 (C2+C3_a0) 


(3.5) 


In  order  to  compute  the  best  two-test  policy  we  use  the  foregoing  one- 
test  solution  for  the  last  step,  and  choose  the  first  test  interval  as  the  x^ 
which  minimizes 


t®/t  \ _ i min 

L2  {T0}  ~ cl  + 


x. 

l 


X 1 

T-1  / ~ [ c?(x.-t)  - c t ] dt 

*0  0 xl  21  3 


+ ( 1 “ Tq  ) [ ” c3  Xj  + Lj  (p  (Tq-  Xj)  ) ] 


(3.6) 


The  bracketed  expression  in  (3.6)  is  a parabolic  function  of  with  a minimum 


at 


c . + ( c - a.) 

*,(2)  * x „ — Tn 

1 C2  + C3"al  0 


(3.7) 
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where 


ai  = P [ a0  + ( c3  - a Q)2/  (c 2 + c3  - aQ)  ] 

However,  if  Tq  < Cj/c^  then  Xj(2)  > Tq  is  not  an  admissible  testing  interval. 

In  such  a case,  the  shape  of  the  x^  - function  indicates  that  6^(2)  = Tq.  It 
follows  that 

62<2>  = T0 

and 

L2°‘T0)  = 2o1  + LS  <V  (3. 8) 

so  the  best  use  of  two  inspections  costs  more  than  using  no  inspections  at  all, 
when  Tq  < c j / c^. 

We  have  just  seen  how  optimization  of  the  first  of  2 tests  might  appear 
to  yield  an  infeasible  solution.  Building  on  this,  a related  difficulty  arises 
when  comparing  different  locations  for  Xj(3)  in  such  cases.  The  remaining 
lifetime  for  the  best  two  tests  may  become  so  small  that  it  is  better  to  re- 
place them  with  no  tests. 

Careful  consideration  of  all  such  details  led  to  an  algorithm  based 
on  interative  calculation  of  the  sequence  of  zero-derivative  intervals  x^(n), 
and  a related  sequence  of  number  Sn.  These  sequences  are  computed  re- 
cursively for  n = 1,2,...  according  to  formulas  given  in  the  appendix,  until 

X1  (»)  > T0  - sn  (3.9) 

* 

If  n is  the  smallest  n to  satisfy  inequality  (3.  9)  then  the  optimal  number  of 
tests  is  found  by  searching  for  the  minimum  of  the  previously  computed  mean 
losses  for  each  number  of  tests  less  than  n . After  thus  determining  the 
optimal  number  of  tests,  their  locations  can  be  computed  from  formulas  given 
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in  the  appendix. 

The  adequacy  of  this  algorithm  is  verified  by  a set  of  lemmas  and 

(5) 

theorems  which  establish: 

(i)  Xj  (n)  > 0 

(ii)  xL  (n)  > Tq  - Sn  implies  L°  > L°_2  (3.10) 

(iii)  Xj  (n)  > Tq  - Sn  implies  x^m)  > TQ  - Sm 

for  all  m > n 


Example 

Figure  3 summarizes  the  res  juts  of  one  example  in  which  the  optimal 
number  of  tests  in  4 (n  =7)  and  the  testing  reduces  the  mean  lifetime  from 
10  to  47.  31. 

4.  Summary 

Tests  which  degrade  performance  have  been  introduced  into  a standard 
model  of  inspection  time  optimization.  Convergent  algorithms  have  been  de- 
veloped for  single-cycle  and  mean  loss  rate  cases  when  the  conditional  life- 
time distributions  are  exponential.  The  uniform  conditional  distribution 
model  was  solved  for  a single-cycle  case.  No  results  have  been  found  con- 
cerning the  rates  of  convergence  of  these  algorithms,  but  all  examples  have 
produced  rapid  convergence,  as  evidenced  by  those  included  here. 

Several  extensions  of  these  results  are  straightforward.  An  algorithm 
for  the  mean  loss  rate  version  of  the  uniform  distribution  case  follows  directly 
from  Theorem  1 and  the  single-cycle  results.  Inclusion  of  a non-zero  pro- 
bability that  the  inspection  causes  immediate  destruction  is  easy.  For 


mm 
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•mmum 


eagsr^  v 


example  the  only  change  thus  introduces  into  (2.6)  is  replacement  of  L.°k+j  by 

< WL°k+1- 

The  general  approach  to  the  exponential  case,  with  backward  recursion 
from  an  approximate  "last-step"  loss  is  applicable  to  other  conditional  life- 
time distribution  of  infinite  support.  Other  distribution  might  require  search- 
type  optimizations  at  each  step,  rather  than  the  simple  ze ro-de rivative  con- 
dition, but  no  other  special  properties  of  the  exponential  distribution  were 
used  in  the  analysis.  This  backward  recursion  might  also  be  efficient  for 
getting  approximate  solutions  to  problems  without  degrading  tests,  by  impos- 
ing a small  mean  lifetime  reduction  factor  in  order  to  have  a "last-step"  loss. 
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7,  Appendix 

This  appendix  lists  formulas  used  to  carry  out  the  algorithm  outlined 
in  Section  3,  all  summations  have  values  of  zero  when  their  lower  limits  ex- 
ceed their  upper  limits. 

a0  = p(c3-C2);  an  = P [an-l+^C3-an-l*  ! (c2+c  3~an-  1 5 n = 1*2*** 
b0  = 0;  bn  (c 3 ” an)  (ncl-iS:ibi)/(c2+C3-an);  n=1’2— 


1q  = °*  = (n  c i "•  ^ p ( c^+c  ^ “ an)  5 n=s  2, 


a mmmmm 


Sn  = [ (n  - 2)  c j - .Sj  b.  + ( c2+c  3-an_2)  S^jj/p  c2 


X1(n)  = [ (n_1)  C1  ' iii  bi+(c3-an-l)Tol  /(c2+c3-an-l) 


n-i-  1 


i “ 1|  2|  • • • n 
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Table  1: 

Single-cycle  Exponential  Case 

Geometric  Failure  Rate  Increase 
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Table  2: 

Single-cycle,  Exponential  Case 

Linear  Failure  Rate  Increase 


t 

c j = 1.0,  = 20.  0,  s = 0. 
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r = 0,  \q  = 5.  0,  M = 21  , p = 0.  9,  k = 20,  and 
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c j = 1.0,  = 20,  s = 0,  r = 0,  p = 0.  9,  M = 2 1,  and  \ = \ q / p* 
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Table  4;  Minimum  Loss  Rates  - Exponential  Case; 


Sensitivity  to  X.q  error 
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and  Optimal  Testing  Times 


